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Abstract 



Using the recently obtained exact results for the expectation values of 
operators in the sine- and sinh-Gordon models [A. B. Zamolodchikov and S. 
Lukyanov, Nucl. Phys. B493, 571 (1997), V. Fateev, S. Lukyanov, A. B. 
Zamolodchikov and Al. B. Zamolodchikov, Phys. Lett. B406, 83 (1997)] 
we calculate the specific heat of the corresponding two dimensional Euclidean 
(classical) models. We show that the temperature dependence of the specific 
heat of the sine-Gordon model, in the commensurate phase, has a maximum 
well below the Kosterlitz-Thouless transition and that the sinh-Gordon model 
is thermodynamically unstable in the strong coupling regime. We give also 
the temperature dependence of the specific heat in the incommensurate phase 
of the sine-Gordon model. 
PACS numbers: 65.40.-f, 65.50. +m 
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I. INTRODUCTION 



The sine-Gordon model is an exactly solvable model which has an enormous number of 
applications in condensed matter physics and statistical mechanics. It has been studied for 
years with many remarkable results being obtained. However, most of the effort has been 
concentrated on studying this model as a quantum field theory. In this paper we discuss 
the sine-Gordon model together with a less famous sinh-Gordon model as models of classical 
statistical mechanics analyzing the behavior of their specific heat in the various range of 
parameters. 

Let us consider the classical sine- and sinh-Gordon models whose partition functions are 
given by 



Z 

'S'sin 
>S'sinh 



Vipe' 



-5M 



T 
T 



sinh I 



(fx 

£x 



|,|V^-Q|^ + ^(l-cos^) 
^|Vv?|V^(cosh(^-r 



(1) 

(2) 



The first model describes, for example, the commensurate-incommensurate transition |T[]. 
Most recently it has been applied to double- layered Quantum Hall systems 0. 

The incommensurate phase appears when |Q| exceeds some critical value and is charac- 
terized by nonzero average value (QVy?). Redefining the field variable {ps/Ty/'^Lp = (f) we 
reduce the above Euclidean action to the canonical sine-Gordon form (see, for example, [^): 
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with = T/ps and /i = m/T, |h| = 27r|Q|//52. 

In a similar fashion the sinh-Gordon action becomes 
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For this model there are no kinks and creation of nonzero field gradient would require imag- 
inary field h. We do not consider such possibility. 



One can consider the quantum field theory in (l+l)-dimensions with Z = j Vcpe'^^'^^ 
where 5* is given by @ or (p. The exact solution for both models is known and we can 
take advantage of the fact that the free energy of a D-dimensional classical model is related 
to the ground state energy Eq of the corresponding quantum field theory living in a space 
of {D — l)-dimensions. More specifically, ioi D = 2 the partition function of the classical 
theory defined on a rectangle x Ly with periodic boundary conditions in the x-direction at 
temperature T is equal to the partition function of the quantum field theory at temperature 
with the coupling constant = T. The limit oo corresponds to the limit of zero 

temperature in the quantum field theory when its free energy is equal to the ground state 
energy Eq = LyE^. Thus we get the following relation between the free energy per unit area 
of the 2-dimensional classical model T and the ground state energy per unit length of the 
(l+l)-dimensional field theory: 



The ground state energy Sq of the quantum sine-Gordon model as a function of parameters 



We start our discussion of the sine-Gordon model with the case Q = which is already 
quite nontrivial. In the Bethe ansatz approach the ground state energy of the quantum 
sine-Gordon model is calculated by regularizing the model by putting it on a lattice. The 
lattice constant a and the inverse coupling constant 9 of the regularized model are related 
to the mass of physical particles. According to 0] the ground state energy (at Q = 0) is 
given by: 



where r = T/{87cps) = T/T^. The parameters 6 and a are related to the kink's mass: 
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/3(T) and is known exactly ( @], 0]) and the corresponding expression for the sinh-Gordon 
model can be extracted from [0. 



II. SINE-GORDON MODEL AT Q = 
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To exclude 6 from Eq. (^) we use the T-dependence of the kink's mass given in 

1 2 r(i^) jnmT{l-^)\^ 



rus = 
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Note that in the hmiting case r ^ 1 we have 

m, = -T-\m/nT,)'/^ . (9) 
a 

We should stress that the above expressions make sense only for msfl <C 1 when the con- 
tinuous approach works. Therefore to calculate the ground state energy in the continuous 
approximation in Eq. (^) one has to keep only the pole closest to real axis. Near the point 
r = 1/2 two poles a.t t = i and t = i/2{l — r) compete and one has to take into account 
both of them. 

Taking this into account we obtain from @ the general expression for the free energy: 

F = F1 + F2 , 
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We emphasize that the necessity to keep both terms in the expression for the free energy 
exists only close to the free fermion point. At r < 1/2 the free energy remains finite in the 
continuous limit (a 0,ms =const), while at r > 1/2 it diverges. The latter fact is in 
agreement with the perturbation theory in m: 

d^x {cos (3(j){x) cos (3(f){0)) ~ J <fx{x/a)-^^^~^^ 

diverges at small distances. 



FIGURES 




FIG. 1. Plots of the specific heat as a function of T/Tc, for different values of the parameter 
m, are shown. 

Since we always keep a finite, we plot the specific heat at finite values of {nm/T^, and it 
is convenient to separate the interval of r into three regions. First, for T/Tc G (0.1, 0.35), we 
take the contribution of only the nearest pole i/2(l— r). Second, the region T/Tf. e (0.35, 0.7) 
where we take the contribution of both poles. Third, for T /Tc G (0.7, 0.9) we take only the 
contribution of the pole t = i. 

We combine these results to find the specific heat C^{T) — —Td^F/d^T as a function of 
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(b) 

m = 0.0000001 




FIG. 2. Plots of {arrLs) as function of (T/Tc), for different values of the parameter m, are shown. 
The condition of continuous approach works on the region {am^) <ti 1. 

T/Tc and m/Tc, in T/Tc G (0, 1). On Fig. 1 we present plots of the temperature of the 
specific heat for various values of m. Fig. 2 gives temperature dependence of the kink's 
mass. From this picture one one can estimate the region where the condition mga <^ 1 is 
fulfilled. 

At 1 — r <^ 1 and at r -C 1 expressions (|l^) for the specific heat simphfy. In the first 
case we have 
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This singularity is associated with the Kosterlitz-Thouless transition at T = Tc 
At r ^ 1 we have 
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The latter expression explains the existence of the maximum in the specific heat: at 
ln(Tc/7rm) ^ 1 the maximum occurs at r* = [ln(Tc/7rm)]~^. 
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III. SINE-GORDON MODEL IN THE INCOMMENSURATE PHASE 



At 

Q>Qc = 4rms(r) (16) 

the sine-Gordon model is in the incommensurate phase characterized by a condensate of kinks 
(QV0) 7^ 0. The ground state energy of the corresponding quantum field theory acquires an 
additional contribution originating from the condensate. The corresponding change in the 
free energy of the classical model is 

SF = P^ + ^ r decos\iet{e) . (17) 

The nonpositive function e(^) is defined inside the interval —B<9<B and satisfies the 
integral equation (see, for example 0) 



t{e) + / de'K{e - 9')e{e') = m, cosh e-— , (is) 

J-B At 
where the Fourier image of the kernel is 



smh L-i { 



2cosh(7rcj/2)sinh2^5^ 

The kernel K{6) encodes the information about the soliton-soliton scattering. The limit B 
is determined by the condition e{±B) = 0. A possibility of e being negative appears when 
the right hand side becomes not positively defined which corresponds to condition (0). 

We have solved the integral equation for the function e{6) numerically and the plots of 
its dependence on r and Q are shown in Fig. |^. The dependence of the limit B on the 
temperature and Q is shown on Fig. |^. 
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FIG. 3. The critical line as function 
of (T/Tc) for a fixed value of the param- 
eter m is shown. The lattice constant a 
is taken to be unity. 

It is curious that the critical field has a finite limit at r ^ (this limit was first considered 
in [^). According to Eq. we have 



aQc(O) = 16 



KnTc 
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In the vicinity of the critical line in the incommensurate phase one can expand solution of 
Eq. (|18D in series in B and get for the additional free energy: 

(20) 
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At small r, KiQ) ^ ^ ln(l/r) and the expansion is valid for 
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Plots of the additional specific heat are shown on Fig. ||a), whereas plots of the total specific 
heat are shown in Fig. ^) for some values of the field Q. 

The Q/t ^ oo analytic structure of the total free energy F{Q) is given by Zamolodchikov 
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The factor k{Q/4T) is given as a power series 
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FIG. 4. (a) Plots of {—e{9)) for different fixed r. On figure (b) the parameter Q = 0.090, wliich 
corresponds to crossing the critical line at r ~ 0. 



Q = ,0600, m = .0001 
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FIG. 5. Plots of temperature dependence of the parameter B in Eq. (| 
taken to have values 0.060 and 0.120, respectively. 



The parameter Q is 
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FIG. 6. (a) Three dimensional plot of the contribution to the specific heat coming from the 
soliton condensate as function of r and Q. The presence of the condensate leads to decrease in the 
specific heat, (b) The total specific heat of the system for different values of the parameter Q. Each 
of the lines 1-5 merges with the line (Q = 0) at the commensurate-incommensurate transition 
temperature. 



The first two coefficients are given by 

, r(r) r(5/2 - r) 
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Putting them together in Eq. ( p2D we get 
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ttKAtJ r(-r)r(l/2 + r) (2r- l)(3-2r) 
The second term on the right hand side of the above equation gives the free energy of the 
system in absence of the field Q (with opposite sign) and cancels with -F(O) of the left hand 
side. On the other hand the ffist and the third terms on the right cancel each other (keeping 
in mind that = Svrpc)- The total free energy of the system in the presence of the field Q 
is given, to this order of approximation, by: 

2 / g\2 r(r) r(5/2-r) r 
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Svr r(-r) r(l/2 + r) (2r - 1)(3 - 2r) 
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The expansion is valid only for small values of r < 1/2, where the free energy of the system 
in absence of the field Q is given by Fi of Eq. (|I0[) . 



IV. SINH-GORDON MODEL; THERMODYNAMIC INSTABILITY 



Now we consider the sinh-Gordon model. Here the exact solution was suggested by 
Fateev who has taken the sine-Gordon two-body S-matrix for the first breathers and 
changed the sign of the coupling constant in it. Comparing Eqs.(^, ^) we see that to get 
the sinh-Gordon action out of the sine-Gordon one we have to change (3 ^ ij3 and m —m. 

Doing this substitution in the expression for the free energy (see also 0), we get 
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where = Svrpg. This expression is valid for T < T^. The integral 
give the expression 
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can be calculated to 
(30) 



with < X < 1. I{x) is monotonically decreasing function taking values in the interval 
(0,1). 

In the light of the following discussion it will be instructive also to have the expression for 
the mass (the inverse correlation length) of the sinh-Gordon theory. To get this expression 
one has to change (3 to ijS for the sine-Gordon mass which corresponds to reversal of the sign 
of Te in Eq.(8): 



Tim r (1 + T/T^ 
~T^T{l-T/T, 



,rfi + ^ T..v(^-TIT^ ' ^^^^ 

(This expression coincides with the mass of the sine-Gordon breather after the substitution 
T — i> — T). We would like to attract the reader's attention to the fact that m(T) never 
diverges at T < and actually goes to zero at T T^,. 
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For r = (1 — T/Tc) ^ 1 we obtain from (p9D /(r) ~ r and substituting this into Eq. p8| ) 
we get at ln(Tc/7rm) ^ 1 

It is easy to see that the specific heat becomes negative at 

(1 - T/Te) < (1 + V2)[2\n{Tj7rm)]-' , (33) 

such that and at temperatures sufficiently close to Tc the model is thermodynamically un- 
stable. 

This thermodynamic instability is not completely unexpected. It occurs in the strong 
coupling regime when the coupling constant T/T^ approaches its critical value 1. Let the 
reader recall that at this value of the coupling constant the ultraviolet limit of the sinh- 
Gordon model - the Liouville model becomes unstable [§ (its central charge approaches the 
value of 25). Another indication of the instability comes from the temperature dependence 
of the inverse correlation length given by Eq. (pT]): at T Tc the mass becomes zero 
which one does not expect to happen to the sinh-Gordon model which can be thought as 
the Gaussian theory perturbed by a strongly relevant operator. One can check that at the 
instability point (|33D the mass of the sinh-Gordon particle is still much larger than its value 
at T = 0. 

V. CONCLUSIONS 

The summary of our results on the classical sine-Gordon model is well represented by 
Fig. 6. The specific heat has a peak well below the Kosterlitz-Thouless transition and the 
temperature dependence becomes even more complicated in the incommensurate phase. We 
suppose that all these features are detectable experimentally in the relevant systems like the 
one described in 0. 
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